Ballistic properties of crystalline defects 
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We introduce a field theoretic formalism enabling the direct study of dislocation and interstitial 
dynamics. Explicit expressions for the energies of such defects are given. We provide links to 
earlier numerical, discrete elastic, time dependent Ginzburg Landau, and other approaches sought 
by numerous authors for the problem of defect dynamics. The formalism introduced in this article 
may be extended to address many other systems. Apart from their heavily studied role in dislocation 
mediated crystal melting, defect dynamics recently gained much interest due to their viable role 
in various electronic and other systems of current interest. These systems include, yet are not 
limited to, stripe phases of Quantum Hall liquids, theories addressing the melting of stripe phases 
in various doped oxides (including the high temperature cuprate superconductors, the nickelates, 
and manganates), defects in supersolids, as well as colloidal crystals. 
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Introduction — The formulation of a field theoretical 
framework for analyzing defect dynamics in two dimen- 
sional crystals is the central goal of this work. Defects 
such as interstitials/ vacancies, dislocations and disclina- 
tions [D, 0, El perturb the ideal perfectly periodic crystal. 
To date, numerous works extensivelystudied the dynam- 
ics of these defects, e.g., 0, 0, 0,11], and some further 
considered how these may, in turn, impede diffusion, e.g. 
Q . Although much progress has been made in analyzing 
defect dynamics, a practical and general field theoreti- 
cal approach that is free from the direct analysis of the 
equations of motion is lacking. Such a general frame- 
work will enable access to numerous physical systems of 
empirical interest. In the current publication, we report 
on precisely such an approach. For simplicity, we focus 
on dislocations and interstitials in two dimensional solids 
free from non-harmonic terms and dissipation. The ex- 
tension to three dimensions is straightforward. The in- 
clusion of disclinations is much more involved. In the two 
dimensional case, the defects are point-like. As these ex- 
citations are well localized, we can wonder if they exhibit 
the behavior of ballistic particles. For instance, what is 
their energy when in motion? The answer to this ques- 
tion is found in Eqs. (f2l| |22|) . our main result showing 
how both interstitials and dislocations run in a 'sound 
barrier' set by the phonon velocities in a way that is 
however different from a simple Lorentz invariance. 

The accumulation of the defects whose dynamics we 
study here is commonly thought to melt the solid phase. 
Early on, Onsager [l0( theorized that the A transition 
in liquid 4 He is driven by an unbinding of vortex lines. 
Shockley [ll[ suggested that the melting of solids, in gen- 
eral, is driven by the proliferation of dislocations. These 
microscopic considerations were couched in a field theo- 
retical formulation by Kleinert Empirical studies of 
these notions have been carried out, e.g., 11211. A scenario 
extending the two dimensional melting of 13] which fur- 



ther suggested an intermediate phase driven by the pro- 
liferation of defects was advanced in 14 , lj| • The glass 



transition has also been suggested to be related to defect 
dynamics H, 17, 18] . In glasses, an extensive configu- 
rational entropy of these defects [l|| (from which ensu- 
ing restrictive slow dynamics might follow (ltSt) m ay be 



sparked. Investigations of Frank Kasper phase s I2C | of de- 



feet lines and of related systems [2JJ, [22J, |23j, |24|, |25 
have flourished. Recent years saw the extension of these 
ideas to various strongly correlated electronic systems 
in which the electronic constituents favored the forma- 
tion of an ideal crystal like stripe pattern. Such patterns 
were detected in the high temperature superconductors 
and other oxides [28[ and in quantum Hall systems [29| . 
Other systems harboring ideal stripe order include fer- 
rofluids, diblock copolymers [3Qj, and colloidal crystals 



311 . 1321 ] . In an interesting development, recent experi- 
ments on solid Helium-4 [33( suggested supersolid-type 
features (a solid analogue of the superfluid phases) 34] 
and have led to a flurry of activity. Elastic defects 3J, 35 j 
may give rise to the observed phenomena either directly 
as superfluid or indirectly [HI]. Indeed, experimental re- 
sults show [37J that the measured putative supersolid- 
type feature is acutely sensitive to the quench rate. The 
results that we report here pertain to defect dynamics in 
these disparate arenas. Specifically, our results hold for 
all classical incarnations or bosonic variants in the quan- 
tum arena of these crystals. In earlier works we employed 
the general framework related to the one used here. This 
enabled us to study other facets of defect dynamics earlier 
and amongst other results, allowed to also find an elastic 
analogue of the Higgs effect 3|| 3t| 4(1 [4l[ and to derive 



the glide constraint of dislocations and to determine new 

H5i- 



non-linear corrections to it 



Elastic gauge fields — Deformations in a two dimen- 
sional crystal are parameterized by a displacement field 
u monitoring the deviations of each constituent "atom" 



2 



relative to its location in the ideal crystal. We begin our 
analysis within Euclidean space (with r = it the imagi- 
nary time) where the Lagrangian density 



(1) 



The Greek indices //, v = 0, 1, 2 with \i = (or v) = 
denoting a temporal component while a,b — 1,2 corre- 
spond to spatial directions alone. To account for the dy- 
namics of the constituent particles, a kinetic term borne 
by C'oo a b = p Sab is inserted (C 0la b = for i ^ 0). The 
spatial components of the elasticity tensor Cij a b host the 
usual harmonic strain energies (V^ — ^diU a Cij a bdjU b ). 
With the kinetic (C TTao ) term in tow, the variational 
equations immediately provide the classical dynamics. In 
the quantum arena, this kinetic term is introduced simi- 
larly for bosonic matter wherein the partition function 



Z = Due 



J d 2 xdrC 



(2) 



We analyze the isotropic case, Cij a b — H>[SijSab + SibSja + 



SinS 



iaVjb\ 



431 ] . At short distances, terms containing 
higher derivatives become significant. For the studied 
isotropic case only two additional terms are symmetry 
allowed. The second order gradient energy density aug- 
menting the first term of Eq. (fTJ) is 



= 2/(V^) 2 + -^£' 2 (V(V ■ u)) 2 



(3) 



The length £, relates to the local rotation (u> = x u) 
inhomogeneities. At distances comparable to length £', 
the compression gradient becomes significant. As cus- 
tomary in studying defects elsewhere, we dualize the ac- 
tion. Stresses are dual to displacements and are given 
by 

Fir 

a ■ ■ r~i o b 

C„ = «TT7^ — = iC^abdvU , 



dC 



ac 



d{d l uj) 



= i4fj£ diui. 



The temporal components <r" are the physical momenta, 
and the spatial erf are elastic stresses. The Hamiltonian 
density, H — ia'Zd® + ir,'9,(V • u) + iridiuo + C, becomes 

n = \ol{c, vah y^i + \\^t>t[ + k^nn. (5) 



A (rotationally related) singularity present in this form 
is readily removed by enforcing the Ehrenfest constraint 



tabO h a 



0, 



(6) 



(with £i2 = —£2i = !)• The dual action is recovered 
from the inverse Legendre transform Cduai = — * cr u^u — 
iT[diiy ■ u) — iTidiUJ + H. The displacement has smooth 
and singular parts u = u smooth + u singu i ar with u singu i ar 
carrying the defect contribution. Upon integration over 



oth 



we are led to the constraint 



(7) 



Eq. ([7]) is none other than an expression of a generalized 
(space-time) conservation law (in the simplest equilib- 
rium setting, the spatial portion of Eq.© dictates a 
divergence-less stress a which implies a net zero force 
acting on each volume element). Eq.© enables us to 
express the brackets as a curl of stress gauge fields, 



-e dr 



(8) 



From Eq.©, it is evident that B is defined only up to a 
gauge transformation, (B" — > BJj + d^f a with arbitrary 
f a leave <r° unchanged). In what follows, we compute in 
the Coulomb gauge diBf — where the calculations sim- 
plify considerably. After the partial integration, the sin- 
gular components become minimally coupled to the dis- 
location currents to the gauge fields, leading to [— i-B" J?J], 
wherein the dislocation currents 



Jl = £ w \d t _ l d\u a . 



(9) 



Following a Fourier transformation, we introduce a new 
basis. For fields of spatial momentum q (and frequency 
bS), we define two planar basis vectors- a longitudinal 
e L = A = (cos 6, sin 6) and a transversal e T = ^4 = 

|q| \ -ri -rj ^ | q | 

( — sin cos (/)) . Expressed in the new basis, the gauge 
fields are B^ = iefBf (double indices are summed over 
with E = L,T) and Bf = -ef ef B^ (with M — r,T). 
In this basis, the Coulomb gauge fix condition becomes 
Bf = 0. Only the longitudinal components of the 
fields t' and r have a physical meaning. This allows 
us to omit the transversal components. Finally, the 
Ehrenfest constraint (Eq.©) may be employed to ex- 
press r = — Bf + i^By, with q = ^\q\ 2 + lj 2 . Armed 
with the above constraints, an integration over the non- 
propagating fields, r' and Bf, leads to 
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+ g 2 (l+f 2 g 2 ) 



1 n^t 

2 p(i + £'V) 

2 T l + £ 2 g 2 



-Bk 



+ q 2 (l+£ 2 q 2 ) 



J 



iB, 



it 



IrH^ 

2 p(l-MV) 

i-i/jJJ+l'VQ/jf- 

1 + f V 

I 



-Bk 



2 r «2 



2/i 1 



iB. 



T\{Jj_ 



1 + f V T 

J L L )+e q \jj + J£) 
l+£ 2 q 2 



(10) 



Eq. (fTU|) constitutes a central result of our work. It en- 
ables a direct evaluation of energy densities from mi- 
croscopic details. We will shortly illustrate how the 
energy densities for different defects (both static and 
dynamic) easily follow from this very general expres- 
sion. Let us first comment on the various terms and 
associated symmetries. The first two terms in Eq. (fT0]) 
contain propagators for the longitudinal and transversal 
"photon" (phonon) respectively. The third and fourth 
terms account for the interaction amongst static dis- 
location charges. Finally, the last two terms embody 
the coupling between the elastic photons and dynami- 
cal dislocation currents. As detailed elsewhere [42j, the 
various current combinations adhere to different rota- 
tional symmetries- (i) the compression (or glide) cur- 
rent [Jj — J^] vanishes for dislocations; (ii) the current 
jj = d T (y x u) — V x (<9 T u) and (iii) the combina- 
tion [jj + Jk} ± i Jt T iJi} i s a snear doublet. We now 
evaluate, by the insertion of the corresponding current 
densities (J) into Eq. (fT0|) . the energies of static and dy- 
namic defects. We start with the static case. Associated 
with a static dislocation of Burgers vector b is a current 
ja;disi. _ 5a<5( r ) Henceforth, we assume that b = be x . 
An interstitial can be viewed as a dislocation dipole - 
a difference of two dislocations carrying Burgers vectors 
of the size of lattice spacing (a) , misplaced by a lattice 
spacing in the perpendicular direction. The current of an 
interstitial J^ mt - = a 2 e a bdbS(r). In Fourier space, 



J 



x\disl. 



J. 



T-int. 



(q) = qa 2 



(11) 



When these external currents are inserted into the action 
of Eq. (fTO)) , we find the static energy of a single disloca- 
tion/interstitial. The range of integration over the mo- 
menta q in all pertinent integrals is 1/L < \q\ < 1/A with 
L the linear dimension of the system size and A the defect 
core size. We facilitate our calculation by setting A = a. 
For a single static interstitial (where an infrared cut-off 
is not imperative), the energy derived from Eq. (|10[) . the 
'action per unit time', is 



T7 1 



27r/ia 



a 2 £' 2 



(12) 



For static dislocations, the energy depends on higher or- 
der couplings and must be regularized by an IR cut-off, 



disl. 




a 2 + £' 2 
1=2. In - 

2 m L 2 



£' 2 



(13) 



The dislocation unbinding energy grows monotonically 
with the increase of the rotational gradient moduli I 
(Eq.©) yet is bounded from above. By contrast, the 
disclination energy Ef^ ic grows with with the square 
of the same coupling. (The latter may be verified by a 
simple rotationally symmetric ansatz.) Thus, high values 
of the rotational stiffness £ suggest a new, intermediate 
phase- in accordance with the known results 14 , 38| . 

We next employ Eq. (fTQ|) to determine the action of 
moving defects. For a uniformly moving defect, the delta 
function constraint of the temporal component of the cur- 
rent is boosted to <5(x — vr) = 5(x — vt)S(u). For the in- 
terstitial motion in all directions is possible. By contrast, 
single dislocation can only_ "glide", i.e. move parallel to 
their Burgers' vectors [l], [J 



38|,|42j. As a consequence, 



in dealing with dislocations we must consider only veloc- 
ities parallel to the Burgers vector. After Fourier trans- 
formation, the temporal currents for an interstitial and a 
dislocation are, respectively, 



J. 



T-int. 



(q) = qa 2 ; J^ sL (q) = M(w - vq x )S(uj - vq x ). 



The corresponding spatial components are 
jT-int. = i vqa ?s(w - V q x ), J%' disl = vbS(uj - vq x ). 

The action per unit time of a kinetic interstitial exceeds 
that of a static interstitial by 



j^int. 1„,2 „„„2 



1 + 



(1 + ^) 2 

2(1+$ 



kin ~ 2 V P na 

2 t n 

2 



(l + z,)-ln(l + - 



(14) 



The analogous 'kinetic' action density for a dislocation is 



rpdisl. 



v 2 pb 2 \ 



111 



L 2 



1 



1 



In 



L 2 



£/2 



(15) 



'Covariant' mapping of a defect gas — With the ac- 
tion densities for both static and dynamic defects at our 
disposal, we may now map our system onto a Ginsburg- 
Landau- Wilson action of the form employed to describe 
a gas of random loops 



= / dxdr 



i|3 M Vf + im 2 H 2 + ^h/f 



(16) 
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With imaginary time as an extra dimension, point de- 
fects become world lines and the action (| 16[) is valid as a 
description of a defect gas. The mapping procedure lead- 
ing to Eq. p^|) , employs a representation of the defect gas 
where, upon conversion with some (yet unknown) veloc- 
ity Cd {xq = CdT), the action of a defect is proportional 
to its world-line length. We will shortly relate Cd to the 
phonon velocity. A defect, moving with velocity v during 
a time At, contributes to the action an amount 



does not lead to an effective dimensionality reduction of 
the defect gas. Instead, a factor of two for velocity in 
Eq. p^]) arises (38[, yet the physical propagators exhibit 
poles corresponding to the physical velocity Cd- 

Similar expressions may be derived within this frame- 
work for single interstitials. The differences are that the 
factor of two lacks from equivalent of Eq. (Tl9]) due to the 
fact that tinterstitials are not constrained by glide, and 
that the interstitial velocity takes form 



/ 2 2 

AS* = e At J 1 + V — « At[ £o + \^—\, (17) 

V C d C d 

with eo the defect core energy. Contrasting this with 
Eqs.CESl US}, we find (for 1,1' < a < L) that 



/' 



I \ 2 i 

\C-d:disl.) — — 



(18) 



Knowing the velocity (Eq. lfTg)) ), the action for the gas of 
dislocations has the derivative term precisely set by the 
first term of Eq. (fT6|) . 

M 2 = Al^| 2 + |V^| 2 . (19) 

The factor of two in the denominator appears for dis- 
location as a consequence of the glide principle- motion 
in a direction transverse to the burgers vector (climb) is 
not allowed for dislocations. Constrained defect motion 



(Cd;i 



H 2(1 + u) 
P l + (2 + vf 



(20) 



Real time ballistics — Upon extending our results to 
real time, we find that the energy diverges as the veloc- 
ity approaches the velocity of 'signal carriers' - phonons, 
v — > ct (the transversal phonon velocity is smaller and 
it bounds the speed of a defect). As the energy has a 
more complicated form than in the imaginary time for- 
mulation, we expose only I = £' = limit where we can 
give exact analytic expressions. With the previously de- 
fined defect currents, now in real time, two gauge fields 
(phonons) will have the solutions 



P 



r 



uJ[), 



p 



-(J T J L 



When these fields are inserted in Eq. (fT0|) and an integra- 
tion over Fourier space is performed, we obtain the total 
(including the static part) energy of defects 



J 



E 



disl. 



(l + z/)(4-3z/) v 2 



2(1-!/) 



1 



, 2l L J 1/(2 - v) 

ir P b 2 In - <^ — T — 

a 1 — v (1 



1 - v 1 - 0- 
I2v{2 - v) 



1 



1 - (3 2 2(1 
6-4/3| 



2 N3/2 



2o^ : 
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(21) 
(22) 



Here, (3t,l = ~^r^- ^s the Lagrangian of Eq.(fT]) is not 
"Lorentz invariant", these expressions do not adhere to 
a single simple covariant form (— 7 ==). Nevertheless, 

we find that divergences occur when /? — > 1. This di- 
vergence of the energy was obtained in our derivation 
by employing the lowest order terms in the elastic gra- 
dient expansion for the Lagrangian density and need not 
constitute a fundamental limitation on interstitial and 
other defect motion. Two velocities appearing in the 
problem complicate the energy further - the low velocity 
approximation yields ballistics (cd) which resemble prop- 
erties from both sectors. On the other side, "causality" 
is driven exclusively by the transverse phonon velocity 
as the energy of a dislocation diverges already at that 



velocity. For small defect velocities, the real time dy- 
namics is not radically different from that in imaginary 
time. At high velocities, retardation effects (triggered 
by finite phonon velocities) come to the fore and the 
(real time) energy vividly displays the role of "causal- 
ity" in the elastic medium. Eq. (|22|) illustrates the power 
of our formalism. With slight variations and different de- 
fect current densities, Eq. (flO|) allows the determination 
of numerous static and dynamic dislocation/interstitial 
configurations. 

Conclusion — In this work, we illustrate that the old 
problem of defect dynamics which is quite cumbersome 
within the common strain formalism becomes far eas- 
ier within the framework introduced here. Much unlike 
the strain formalism, where higher order couplings are 
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necessary to close a set of differential equations, in the 
formalism studied in the current article, we obtain ex- 
pressions which are exact and fairly trivially. To illus- 
trate the power of our formulation, we compute the en- 
ergies of moving dislocations and interstitials (Eq.(|22p). 
With slight variations and different realizations of defect 
current densities, Eq. (|10p allows the analysis of numer- 
ous static and dynamic dislocation/interstitial configu- 
rations. As the problem of elastic defect dynamics has 
myriad manifestations - crystals, glasses, liquid crystals, 
several strongly correlated electronic systems, and many 
others - our study affords a fresh perspective on the de- 
fect dynamics within these systems. 

We conclude with two brief remarks — (i) The dynam- 
ics found in imaginary time has a crucial consequence on 
various properties of a special nematic phase of a crys- 
tal wherein only dislocations proliferate [lij]. Here, dy- 
namic dislocations alter physical measurable quantities 
such as elastic response as well as electric response func- 
tions (4oj ]. (ii) Within this work, only dislocation defect 
currents were treated. To accommodate disclination de- 
fects, we need to define the fields B of Eq.® such that 
B l = 9 a A^ + e ab d b Al. The "double curl gauge fields" 
Aj;' T are defined up to innocuous shifts, 8A^' T — d^g L ' T . 
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